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Abstract 

At future electron-positron colliders, one of the largest irreducible backgrounds to top 
searches in the channel '4 jets + lepton + missing energy' comes from QCD events of 
order a 2 s . We compute here such processes exactly at the parton level by resorting to 
2^6 matrix elements exploiting helicity amplitude techniques. We adopt a typical 
selection procedure based on the tagging of a high momentum and separated lepton. 
We finally outline kinematic differences between signal and background events that 
can be exploited to further reduce such a QCD noise. 
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1 Introduction 



One of the top on the list reasons to build an e + e~ linear collider (NLC) with a centre-of- 
mass (CM) energy between 350 and 1 TeV [|T| is to study in great detail the top parameters: 
its mass (m t ), width (T t ), quantum numbers (Q t , If) and branching ratios (BRs) 0. Not 
surprisingly so, as it is not unreasonable to believe that the heaviest of the fundamental 
particles discovered so far would after all have something to teach us ||. Indeed, in the 
unforeseen scenario that no Higgs bosons and no Supersymmetric (SUSY) particles are 
found at the LHC [|J, this might even be the only task left for the NLC: to run as a top 
factory^]. Needless to say, under such a gloomy prospect, one would want to make the most 
out of such a machine, one way or another. 

Our contribution in that respect is to calculate the 2 — > 6 scattering processes 

e + e- -> qq QQ' £u e (1) 

e + e - ^ gg QQ' i Uf (2) 

at the parton level, through the perturbative order 0(a^ m a^), for any possible flavour 
combination of quarks q,Q,Q' = u,...b ^ t and leptons £, z/£, with £ — e, \i ^ r. The 
sum of these two mechanisms represents one of the largest 'irreducible' backgrounds to top 
production and decay in the semileptonic channel 

e+e" -^ti^bb QQ' Zv h (3) 

the one preferred for experimental studies ||. Their calculation has never been attempted 
before^. Obviously, if one wants to perform precision measurements of top parameters at the 
NLC, to pin down the size and shape of all important background processes is of paramount 
importance. 

The plan of the paper is as follows. Next Section describes our calculations. The third 
one presents our results. Conclusions are in Sect. f|. 

2 Calculation 

Most of the Feynman diagrams that one has to tackle in order to calculate processes (0)-(H) 
proceed through + 4 jet production |J, with the gauge bosons subsequently decaying to 
lepton-neutrino pairs, see first four (six) graphs in Fig. |l](@). In addition to these, one also has 
to consider the graphs in which the off-shell W^* boson is produced in the bremsstrahlung 
off a leptonic current, see the last one(two) in Fig. |l|(@), and eventually decays to four 
partons. Altogether one has to compute 140 diagrams: 32 associated with process (|1|) and 
108 with (fj). (Notice that for the time being we neglect diagrams in which an electron 

1 The other, more frightening legacy of a Higgs- and SUSY- less LHC would be a strongly interactive weak 
sector, building up somewhere around the TeV scale, the very upper end of, if not beyond, the technical 
reach of the NLC as well as of many perturbative calculations ! 

2 A preliminary exercise in such direction was performed in |^|, where however only the case of on-shell 
production plus four jets was considered. For the case of 0(a® m ) backgrounds see Ref. 0. 
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and the companion neutrino in the final states of ([I])-(|2]), i.e., i = e, are connected to the 
incoming beams via so-called 'multi-peripheral' channels. We will come back to this point 
in the conclusive Section.) The signal comes via two simple s-channel graphs (that we do 
not reproduce here). 

We calculate the signal (|3|) at the leading-order (LO), though we are aware that sev- 
eral higher order electroweak and QCD corrections (mainly to the on-shell production) are 
known to date || f|. We do this for consistency, as the background processes ©-(1) can 
only be evaluated at tree-level with present technology. To compute all graphs is not a 
prohibitive task, if one resorts to helicity amplitudes methods. We have done so, by using 
both MadGraph || and a self-made program based on the technique of Ref. [|1(|. They 



agree with each other. Moreover, they have passed all our gauge-invariance tests, so to give 
us confidence in our numerical results. 

To obtain the latter, one has to integrate the Feynman amplitudes squared over a six- 
body phase space. This task is not difficult either, provided one takes some special care 
in dealing with the various resonances. To get around this problem we have proceeded as 
described in Ref. [IT]. That is, by splitting the gauge- invariant matrix elements of ([]])- 



(0) in several sub-terms, each having its peculiar resonant structure. To any of these a 
dedicated mapping of the phase space has been attributed. These pieces have eventually 
been summed up together, after integration (performed with different packages, for cross- 
check purposes), so to recover gauge-invariance. A flat phase space has been instead used 
to integrate the interferences between the different sub-terms, with the help of some brute 
force too (i.e., a largely increased number of random calls). In general, we have verified that 
their contribution is never dominant, but not necessarily negligible, as compared to the pure 
resonances. For reason of space, we will not dwell here in technicalities any further, as we 
will discuss only the total integrated rates, summed over all production sub-channels. 

Before proceeding to present our results, we list the numerical values adopted for the 
various Standard Model parameters: 

mi = m n = m u = m d = m s = m c = 0, 

m b = 4.25 GeV, m t = 175 GeV, 

M z = 91.175 GeV, T z = 2.5 GeV, 

M w = 80.23 GeV, T w = 2.08 GeV. 

As for the top width T t , we have used the leading-order (LO) value of 1.55 GeV as a default. 
However, in a few cases, we have compared the yield of process (Q) with finite width effects 
to that of the same reaction in Narrow Width Approximation (NWA), for which we have 
rewritten the (denominator of the) top quark propagator as 

i / r \ 1/2 

(4) 



p 2 — m 2 + im t T \T t 



with r — > 0, thus mimicking a delta function?] 5(p L — m 

3 For r ~ 10 -5 GeV the total cross sections in NWA are stable and reproduce the on-shell results of 
the 2^2 process e + e~ — > tt within numerical accuracy. (Note that for r = Tt in eq. (Q) the standard 
expression of the propagator is recovered.) 



2 



For the vector and axial couplings of the gauge bosons to the fermions, we use the 
'effective leptonic' value 

sm 2 (6 w ) = sm 2 e{i (6 w ) = 0.2320. (5) 

The strong coupling constant a s entering processes (fl])-(|2]) has been evaluated at two loops, 
with Nf = 5 and A^jg = 200 MeV, at the scale Q 2 = s. The electromagnetic one was set at 
1/128. Finally, the centre-of-mass (CM) energies considered for the NLC are yfs = E cm = 
360 and 500 GeV, as representative of the threshold, y/s ~ 2m t , and asymptotic, y/s ^> 2m t , 
top-antitop production regimes. Beyond those energies is no longer the top quark realm. 

3 Results 

In discussing the interplay between the background (|Tj)-(F2|) and signal (|^) processes, we 
have focused our attention to the case of the semi-leptonic (or, equivalently, semi-hadronic) 
signature 

4 jets + (6) 

where £ = e or /x, ]p represents the missing energy/momentum due to the neutrino vg 
escaping detection and where to the four-jet hadronic system no 6-tagging is applied^. We 
identify the quarks and gluons in with the jets in @ and apply all our cuts directly 

at the partonic level. 

The choice of considering here only the final state (|6]) is indeed not restrictive, in the sense 
that the latter is to date the experimentally preferred channel in searching for e + e~ — > ti — * 
bbW + W~ events. On the one hand, as opposed to the fully hadronic signature bbW + W~ — > 
6 jets, it has a simpler detector topology and thus is much easier to reconstruct, further 
allowing one to reduce the severe problems due to the six-jet combinatorics (particularly, 
if no heavy quark identification is exploited)^. On the other hand, the case involving two 
leptonic W ± decays has a double disadvantage as compared to channel ([]), that is, a very 
much reduced statistics and problems in reconstructing the top mass spectra because of 
the two neutrinos. Furthermore, as selection method of candidate top-antitop events we 
adopt one rather similar to that outlined in Sect. 4.2 of Ref. ||, based on the detection of 
high-momentum isolated leptons. As a matter of fact, such a procedure has been shown 
to be the most effective one, as it eventually yields the largest signal-to-background ratio, 
both at and above threshold: see Tab. 4.2 of |J. 

About 40% of top-antitop events (^) produce an energetic electron or muon which is 
clearly separated from the hadronic system. Following Ref. ||, we consider a lepton to be 
isolated if a jet clustering algorithm with a 'minimum mass' cut-off recognises it as a 'jet' 
with a single particle. As jet finder we use the Jade one [fL^| , with y > 0.003(0.002) at 
yfe = 360(500) GeV, so that M £±j = 2E i ±E j (l - cos6 l±j ) > 19(22) GeV for each jet j in 

4 If the latter is enforced, we have already shown in Ref. |(| that the irreducible background to top events 
due to l W ± + 26 + 2 jet' events is negligible, provided that a vertex tagging efficiency ~ 0.5 can be 
achieved. 

5 Also notice that, once the missing energy /momentum has been assigned to the neutrino, the kinematics 
of (p|) is fully constrained, like in the case of the purely hadronic channel. 
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(P). As the MEs of the background processes are divergent if the partons are allowed to be 
infinitely soft and/or collinear, we also apply the jet clustering algorithm to the hadronic 
part of the event, for all cases ([]])-©. Fig. | presents the y-dependent total rates for these 
three processes as obtained by enforcing the jet clustering algorithm only, e.g., at = 360 
GeV. Notice that there exists a hierarchy in the production rates: 

a(e + e~ ->• it -> bb QQ' Ivi) > (x(e+e~ -> gg QQ' £u e ) > a(e + e~ -> qq QQ' £v e ). (7) 

At the minimum value of y considered here, they approximately scale as 100:10:1. If y/s = 
500 GeV, see Fig. |], the relative ratio of process ([|) to either of (|2|) or (Q) further increases, 
while that between the latter two suffers little from the CM energy scaling. As for top width 
effects, whereas these are naturally sizable at threshold (with differences of about 10%, see 
the top-right insert in Fig. |]), they instead fall to the percent level in the asymptotic regime 
(Fig. § JTJ 



In Figs. we have neglected considering Initial State Radiation (ISR) fllijl , that is, 
the presence of bremsstrahlung photons emitted by the incoming electron/positron beamsQ. 
The main consequence of ISR is to lower the effective CM energy, thus ultimately reduc- 
ing(enhancing) the production rates of processes whose cross sections increase (decrease) 
with yfs. One thus expects the top-antitop rates (|3|) to be rather sensitive to ISR, for 
two reasons. Firstly, because of the s-channel topology of the Feynman diagrams involved 
(which tends to increase the rates). Secondly, at threshold, because the difference \fs — 2m t 
starts approaching the edge of the phase space (thus decreasing the rates). In contrast, the 
background rates ([I])-(|2]) should depend much less on the ISR. On the one hand, they are 
not purely s-channel. On the other hand, \fs is well above the heavy particle thresholds 
which can onset there (such as the dominant W 7 " 1 "!^ - ). 

This dynamics can be appreciated by comparing the total cross sections in the upper lines 
of Tab. |l|, when no selection cuts are applied apart from the jet clustering algorithm. There, 
notice that the background rates (third and fourth column) outside and inside brackets are 
rather steady, at both collider energies. For signal events, both in NWA and with finite 
top width (first and second column), differences are much more sizable. At threshold, it is 
clearly the phase space suppression to dominate, depleting the signal rates by up to 43% (at 
the minimum y). In the asymptotic regime, where phase space effects become negligible, 
the s-channel increase is overturned by the presence of the invariant mass constraints, as 
also typical energies of the final state particles diminish because of ISR. 

Indeed, in presence of the latter, the response of processes (0)-(@) to the implementation 
of any selection cuts is no longer straightforward, as a consequence of the fact that ISR also 
induces a smearing of the differential distributions. Thus, from now onwards, all our results 
will include initial state bremsstrahlung. (We will keep those without it only for reference 
purposes.) Among the various ways of implementing the ISR, we have adopted here the so- 



called Electron Structure Function (ESF) approach, based on the formulae given in Ref. [[15 
In addition, hereafter, we will stop considering the case of process (|3]) in NWA. 



6 Also beamsstrahlung and Linac energy spread E3 in principle affect processes In practise, for 

narrow beam designs of the NLC, their effects are much smaller as compared to those induced by ISR H], 
so for the time being we neglect the former in our analysis. 
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We now proceed by applying all other selection cuts of || that can be exploited at parton 
level too. Namely, an event is accepted if: 

1. its 'thrust' (calculated by using the four jet and the lepton momenta), T e ± t j, is sig- 
nificantly far from the infrared region typical of QCD events; 

2. the invariant mass of the hadronic system, M 4j -, is far above the typical resonances of 
background events (Mw in our case); 

3. the amount of missing energy, £*, is rather contained, as in top events this is typically 
less than \jm\ — M^; 

4. the (absolute) momentum of the isolated lepton, \p£±\, is above a minimum energy 
threshold and below a maximum one of standard acceptance. 

Numerically, to account for other sources of background too, other than ([1])-(|2]), we require 

T e ± t j < 0.75 M 4j > 0.4v^ 
$ < 0.4^ 0.04^ < \p e ± \ < 0.3y/s. (8) 

The second line in Tab. [l] reports our findings. After the cuts in (H) are enforced, the 
background from processes (|I|)-(fJ) amounts to about 1.6% of the signal @ at threshold, 
whereas well above that the corresponding figure is ~ 0.9%. Thus, in both cases, the QCD 
noise is under control. 

Nonetheless, one ought to know its effects on the differential spectra used to fit the top 
parameters. We consider here all possible three-jet mass distributions which can be 
reconstructed in samples of the type (|6|). After ordering the four jets in energy, such that 
E\ > ... > £4, one can build up four ijk combinations, such that i < j < k. They are 
presented in Figs. || and ||, for the cases y/s = 360 and 500 GeV, respectively. The size and 
shape of the backgrounds ([l])-(|2]) are rather innocuous in the vicinity of the top peaks, so 
that one should not expect any significant distortion of the Breit-Wigner distribution of the 
top resonances. As a matter of fact, in this respect, it is the intrinsic background due to 
mis-assigned jets originating in the signal ([3]) from 6-quarks that affects most the signal, as 



discussed in Refs. 11, 13 



Finally notice that an additional requirement can be imposed to events of the form 
(H), in order to increase the signal-to-background ratio of @ vs. (|l])-(@). That is, that 
one two-jet combination ij, among the six possible possible ones, when i < j = 1, ...4, 
produces an invariant mass around the W ± mass. As one can appreciate in Figs. |7|-[8], 
this would reduce the QCD noise to imperceptible levels. If one imposes at ^/s = 360 
GeV, e.g., |M 34 — M w \ < 20 GeV, than additional reduction factors of 9.1 and 8.4 apply 
to the processes ([!]) and @, respectively, whereas the loss on the signal (0) is just 1.3. 
Corresponding numbers at yfs = 500 GeV are 4.3 and 3.7 for the backgrounds, and 1.6 for 
the signal. Finally, we should mention that we have tried other quantities too (such as, e.g., 
jet energies, relative angles, etc.) but they have proved themselves much less useful than 
the Mij spectra in disentangling reactions (|I])-(Q) and @. 
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4 Conclusions 



Thus we conclude that irreducible 0{ot\ m a 2 s ) backgrounds to the '4 jets + lepton + missing 
energy' signature of top-antitop events at the Next Linear Collider are reduced at the ten 
percent level by using a standard selection procedure (at typical design energies), in line 
with previous results obtained by using parton shower models. Such figure can vigorously 
be reduced further if a simple mass requirement on a two-jet system is imposed. We have 
obtained such results by computing tree-level matrix elements at leading-order for the rel- 
evant 2^6 processes, with the only exception of multi-peripheral graphs entering final 
states including electrons. We have neglected the letter for two reasons. On the one hand, 
we would have had to calculate twice as many diagrams as compared to the way we did it. 
On the other hand, the contributions of the missing terms has already been proved to be 
very small in the case of 0(a & em ) e + e~ — > bb QQ' £u£ processes, as their inclusion account 
for an increase of only 6% with respect to the muon rates, for a CM energy of 500 GeV 
[0]. Indeed, we expect the same to have occurred here. Anyhow, given the results we have 
eventually obtained for the signal-to-background ratio, even if their rate is actually much 
larger than the mentioned figure, our conclusions would remain unchanged. Finally, we be- 
lieve that, although confined at the partonic stage, our findings should not be invalidated by 
studies at the hadronic level. We make our programs available to the public for simulation 
purposes in the above respect. 
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e+e" -> 4 jets + £± + $ at the NLC 


o"tot (fb) 


ft (NWA) 


tt -> bbw+w 






v/i = 360 GeV, % J ± . > 0.003 (My > 19 GeV) 


76(53) 
40(28) 


70(49) 
37(26) 


6.5(6.4) 
0.50(0.38) 


0.42(0.42) 
0.024(0.023) 


yfi = 500 GeV, % J ± • > 0.002 (My > 22 GeV) 


102(99) 
32(36) 


102(99) 
33(36) 


3.7(3.7) 
0.55(0.32) 


0.21(0.20) 
0.025(0.019) 


T <±i j < 0.75 M 4i 
^ < 0.4^/i QM^/s <\p e ±\ 


> 0.4^ 
< 0.3v/i 


No 6- quark tagging 



Table 1: Cross sections of processes fl]])-®, the latter in both NWA and with finite width, 
at y/s = 360 and 500 GeV. First line is without the kinematical cuts (|8|). Second one is with 
the latter implemented. In brackets, the same rates in presence of Initial State Radiation. 
A default jet clustering procedure is enforced in all cases. 
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Figure 1: Relevant Feynman sub-diagrams contributing at lowest order to process ([!]). 
Permutations of real and virtual lines along the fermion lines are not shown. An internal 
wavy line represents a W , a 7 or a Z, as appropriate. 



9 




Figure 2: Relevant Feynman sub-diagrams contributing at lowest order to process (@). 
Permutations of real and virtual lines along the fermion lines are not shown. An internal 
wavy line represents a W , a 7 or a Z, as appropriate. 
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Figure 3: Cross sections of the processes: (|3|) in NWA (solid) and with finite top width 
(dashed), (|) (dotted) and (0) (dot-dashed), function of y for the Jade algorithm 

(applied to the four jets and the lepton), at y/s = 360 GeV, before the selection cuts (§) 
and without ISR. The insert refers to the top-antitop rates only (labelled as above). 
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e e -» 4 jets 1 E(miss) 
E cm = 500 GeV ' No cuts 

1 f i i i i i i i i I i i i i ; 




0.005 0.01 0.05 
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Figure 4: Cross sections of the processes: (|3|) in both NWA and with finite top width (solid, 
the two visually coincide), @ (dotted) and (jl]) (dot-dashed), as a function of y for the Jade 
algorithm (applied to the four jets and the lepton), at y/s = 500 GeV, before the selection 
cuts (H) and without ISR. 
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e e -» 4 jets 1 E(miss) 
E cm = 360 GeV, y(J) = 0.003, cuts 
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Figure 5: Differential distributions in the invariant mass of the energy-ordered three-jet 
pairs, M^k with i < j < k = 1, ...4, for the processes: (|3]) with finite top width (solid), 
(H) (dotted) and (0) (dashed), for y = 0.003 in the Jade algorithm (applied to the four jets 
and the lepton), at yfs = 360 GeV. The selection cuts (|8|) have been enforced and the ISR 
implemented. 
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e e -» 4 jets 1 E(miss) 
E cm = 500 GeV, y(J) = 0.002, cuts 




M ijk (GeV) 



Figure 6: Differential distributions in the invariant mass of the energy-ordered three-jet 
pairs, M^k with i < j < k = 1,...4, for the processes: (|3]) with finite top width (solid), 
(H) (dotted) and (U) (dashed), for y = 0.002 in the Jade algorithm (applied to the four jets 
and the lepton), at yfs = 500 GeV. The selection cuts @ have been enforced and the ISR 
implemented. 
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e e -> 4 jets 1 E(miss) 
E om = 360 GeV, y(J) = 0.003, cuts 
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Figure 7: Differential distributions in the invariant mass of the energy-ordered two-jet pairs, 
My with i < j — 1, ...4, for the processes: with finite top width (solid), (Q) (dotted) and 
(P (dashed), for y = 0.003 in the Jade algorithm (applied to the four jets and the lepton), 
at y/s = 360 GeV. The selection cuts (§) have been enforced and the ISR implemented. 
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Figure 8: Differential distributions in the invariant mass of the energy-ordered two-jet pairs, 
My with i < j — 1, ...4, for the processes: with finite top width (solid), (Q) (dotted) and 
(P (dashed), for y = 0.002 in the Jade algorithm (applied to the four jets and the lepton), 
at y/s = 500 GeV. The selection cuts (§) have been enforced and the ISR implemented. 



16 



